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ABSTRACT 


The  relationship  between  the  stress  relaxation  modulus,  E(t) , 

dynamic  modulus,  E*(icu)  ==  E*  (cu)  +  i  E”(cu),  can  be  expressed  in  the 

00 


E*(a)) 


and 


E"(£u) 


)  =a;  fl 

Jo 

=  w  /E 

Jo 


E(t)  sin  ojt  dt 


00 


E(t)  cos  tut  dt 


and  the 
form 


The  evaluation  of  these  integrals  must  be  a  numerical  process  since  E(t) 
as  obtained  from  test  data  is  ordinarily  given  in  tabular  form.  The  basic 
numerical  approach  presented  in  this  document  is  to  fit  E(t)  over  each 
time  definition  interval  with  a  function  of  the  type 


E(t) 


'n-Hl 


=  A  e 
n 


-b  t 
n 


t 

r 


The  integral  defining  E*(cu)  can  then  be  expressed  as  sums  of  integrals  of 
the  type 


^  n+1 


J 


t 

n 


A  e 
n 


sin  cut  dt 


This  integral  and  the  companion  integral  for  E"(iu)  can  be  evaluated  analyt¬ 
ically,  resulting  in  a  series  approximation  for  E*(i<u)  containing  no 
integrals . 


The  numerical  procedure  is  carried  out  using  a  digital  computer.  An 
accuracy  investigation  is  performed  using  the  tensile  stress  relaxation 
data  for  polyisobutelene.  The  regime  in  which  E(t)  must  be  known  in  order 
to  accurately  compute  E*(icu)  is  determined  for  a  frequency  spectrum. 
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SECTION  I 


INTRODUCTION 

The  stress  relaxation  modulus,  ECt),  for  a  linear  viscoelastic  mate¬ 
rial  can  be  converted  to  the  dynamic  modulus,  E*  (i(u)=E '  (cu)+iE"(aj) ,  by 
evaluating  the  following  integrals:^ 


/.OO 

E '  (w)  =  CJ 

/  E(t)  sin 

Jo 

cut  dt 

^00 

E'(aj)  = 

1  E(t)  cos 

Jo 

cut  dt 

It  is  not  possible  to  obtain  E(t)  from  tests  in  the  time  range  speci¬ 
fied  by  either  the  lower  or  upper  limit  of  these  integrals.  Therefore, 
what  must  the  time  range  be  for  E(t)  in  order  to  obtain  a  good  approxi¬ 
mation  for  E*(ia;)?  The  data  from  stress  relaxation  tests  is  ordinarily 
obtained  in  tabular  form.  A  numerical  procedure  must  be  used  to  evaluate 
the  integrals  for  E*  (cj)  and  E”(cu)  . 
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SECTION  II 


METHOD  OF  SOLUTION 

Given  the  time  dependent  relaxation  modulus,  E(t),  for  a  linear 
viscoelastic  material,  it  is  desired  to  compute  the  frequency  dependent 


dynamic  modulus,  E*(ia;)  =  E’(aj) 

+  i  E"(co). 

The  relationship  between 

E*(iu))  and  E(t)  can  be  expressed 

in  several 

forms.  Commonly  used  forms 

for  this  relationship  are : 1 

^  00 

E '  (cu)  =*  CO 

1  E(t)  sin 
•'o 

wt  dt 

(1) 

^  00 

E"((o)  =»  (0 

1  E(t)  cos 

Jo 

ojt  dt 

(2) 

and 

/•OO 

E  *  (co)  »  E^  4-  CO 

/  \/^(t)  sin 

wt  dt 

(3) 

•'o 

^00 

E"(co)  =»  CO 

/  \/^(t)  cos 

cut  dt 

(4) 

• 

^0 

where:  '//'(i)  =  E(t)  - 

=  constant  value  approached  by  E(t) 

It  can  be  shown  that  equations  (1)  and  (3)  are  equivalent,  as  are  equations 
(2)  and  (4)  .  This  demonstration  is  based  upon  the  fact  that 

^00  ^00 

/  sin  X  dx  =  1  and  |  cos  x  dx  =  0 

Jo  Jo 

The  integral  in  equation  (1)  can  be  written  as 


sin  cjt  dt 


2m  TT 

(0 

.00 

E(t)  sin  (ot  dt  + 

1  E  sin  cot  dt 

J  ^ 

2tn_7r 

where  m  is  an  integer  chosen  such  that  E(t) 


E  for  all  t>  2m 7r 
^  0) 


1,2 
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Since 


/■ 


00 


u)  I  sin  cut  dt  => 

ZltlTT 

cu 


and 


Zmir  .  TT 
cj  2  cu 

sin  cut  dt 
c 

2in7r 

cu 


equation  (1)  can  be  written  as 


E'(a,) 


'  2m  TT  +  TT 
2  cu 

^  cu  I  ^(t)  sin  cut  dt . 


Z.IU  7 

t  " 

/  E(t 
•'0 


A  similar  treatment  of  equation  (2)  yields 


2m  TT 
a) 


E*'(cu)  =  cu  /  E(t)  cos  cut  dt, 

■^0 


(5) 


(6) 


Because  E(t),  as  obtained  from  tests,  is  in  tabular  form,  a  numerical 
method  must  be  used  to  evaluate  the  integrals  for  E*(ia^ .  A  procedure  for 
curve  fitting  will  be  necessary  in  any  method  chosen.  Two  distinct  types 
of  curve  fitting  methods  are  available.  One  type  consists  of  determining 
a  single  function  which  is  good  for  the  full  range  of  the  independent 
variable.  This  is  the  type  of  fit  accomplished  when  a  multiparameter 
Maxwell  model  is  assumed  for  the  material.^  The  other  type  of  curve  fit 
method  consists  of  selecting  discrete  functions  to  be  used  in  different 
ranges  of  the  independent  variable.  The  latter  method  was  chosen  for 
further  study  because  of  its  simplicity  and  accuracy  potential. 

Examination  of  several  sets  of  tabular  data  for  E(t)  has  indicated 
that  E(t)  can  be  fit  over  time  intervals  with  functions  of  the  type 


E(t) 


n-fl 
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(as  shown  in  Figure  1) 


Interval  Between  10  and  100  Hours 


Substitution  of  this  expression  into  equations  (5)  and  (6)  yields 


E'  (oj)  =  o) 


“1 

I 


“n+l 


-b  t 
n 


sin  cjt  dt 


n=l  n 


(7) 


E”(cj)  =  0) 


n+1 


cos  ojt  dt 


(8) 


n=l  n 


where  is  chosen  such  that 


and  n2  is  chosen  such  that 


2in7r  jr^ 

u)  2^ 


2njjr 


For  the  time  interval  between  t^  and  integrals  can  be  evaluated 

analytically,  the  result  being 


JL 

\  an 

E '  (to)  =0)  )  ,  2+  2 

/  .  n 


n-1 


^  t_ 


e  (b^  sin  oit  +  oi  cos  wt) 


n 


n+1 


(9) 


and 


E**(a))  =*0) 


n 


u  2 

b  +UJ 
n 


n=l 


e  (b  cos  CJt  -  01  sin  cjt) 

^  n 


(10) 


-I  ^n+1 


I 


1 
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where: 


w  =  given  circular  frequency 


b  =  In 
n 


a 


n 


n+1 

E(tn) 


t 

n 


Two  considerations  should  be  kept  in  mind  when  the  time  increments,  t  ,  in 
equations  (9)  and  (10)  are  chosen .  The  first  consideration  is  that  for 
small  values  of  t,  E(t)  changes  rapidly  with  time.  The  t^  should  be 
chosen  in  this  time  regime  where  the  test  data  points  occur.  The  second 
consideration  concerns  convergence . 

A  method  of  determining  the  convergence  of  the  integrals  in  equations 
(5)  and  (6)  should  be  provided.  A  convergence  criterion  is  easily  obtained 
if  t^+i  is  chosen  such  that  sin  cut^+l  =  sin  cl>  ^n  and  cos  =  cos  . 

If  tn  is  large  compared  with  27r/cj  ,  then  the  data  points  thereafter  can  be 
assumed  to  occur  at  multiples  of  27r/aj  with  little  error. 

Assume  t^-fl  =  t^  +  2  ^  in  equation  (9)  j  where  p  is  a  given  positive 
integer.  Further,  assume  that  a^^,  b^,  and  oi  are  given  for  the  interval 
tn  to  Then  the  interval  contribution,  8^,  is  equal  to 

t 

I  n 


t  2p7r 
^  CO 


Sn 


~h 


2^ 

U) 


b  2-t-a;2 


“bnt 


^  (b^  sin  ojt  +  cos  <i>t) 


This  equation  can  be  expressed  in  the  form 

221 


(1-e 


-Ipw 


^ )  I 


sin  (a)  tn  +  ^  )[  (11) 


where  ~  tan 


-1 


U) 

(r“> 


( 


I 
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A  plot  of  83^  as  function  of  t^^  will  appear  as 


TT  ^ 

The  interval  contribution  is  a  maximum  when  t-  =  q  t —  -  ,  when 

^  2  cj  ^ 

q  IS  an  even  integer j  and  the  contribution  is  zero  when  q  is  an  odd 
integer . 

As  E(t)  approaches  a  constant,  bj^  and  cf>  approach  zero*  Thus  8^  is 
small  when  t^  is  an  odd  multiple  of  wlloj.  This  fact,  coupled  with  the 
fact  that  t^^  must  be  an  odd  multiple  of  ,  has  the  implication  that 

the  t^'  in  equation  (9)  should  be  chosen  as  odd  multiples  of  7r/2cj  when  t^ 
is  sufficiently  large  to  permit  the  approximation.  With  the  t^  thus 
chosen,  equation  (11)  reduces  to 


A  similar  discussion  and  treatment  of  E"(cj)  implies  that  t^  in  equa¬ 
tion  (10)  should  be  chosen  as  even  multiples  of  7r/2cu  when  t^  is  suffi¬ 
ciently  large  to  permit  the  approximation.  With  these  assumptions 
concerning  t^^,  the  interval  contribution  formula  for  E”(cu)  is  identical 
with  equation  (12). 

A  computer  program  was  written  for  the  IBM  1620  digital  computer. 

(See  Appendix  A.)  This  program  evaluates  equations  (9),  (10),  and  (12). 
The  time  intervals,  t^  and  In  equations  (9)  and  (10)  are  chosen  as 

the  intervals  at  which  the  data  points  are  given.  When  a  t^  is  found  such 


7 


that  K  ^ — jjj — j  5  where  K  is  an  integer  given  as  input,  the  program 

assumes  that  E(t  )  is  equal  to  E(t  )  where  ojt  =  (  ^  +  nearest  integral 
n  P  P 

multiple  of  Itt)  in  equation  (9)  and  ut^  =  nearest  integral  multiple  of 
2 TT  in  equation  (10).  Thereafter,  the  program  uses  equation  12  to  compute 
E*(icu)  . 
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SECTION  III 
DISCUSSION  OF  RESULTS 


To  analyze  the  accuracy  of  the  method,  the  stress  relaxation  data  for 

4 

polyisobutelene  (see  table  I)  was  run  using  this  computer  program.  The 

stress  relaxation  data  extends  over  16  decades  in  time,  and  is  acknowledged 

to  be  the  best  data  available  at  the  present  time.  The  computer  results 

for  E*  (ioi)  are  shown  in  figure  2.  The  data  for  the  dynamic  shear  modulus, 

5 

(icu) ,  for  polyisobutelene  was  converted  to  E*  (icu)  using  the  formula 
E*  (ioj)  =  3  G*  (icu).  The  results  are  also  plotted  in  figure  2.  The  com¬ 
parison  is  excellent  throughout  most  of  the  frequency  regime. 

The  data  was  further  analyzed  to  determine  the  sum  interval  contri¬ 
butions  in  equations  (9),  (10)  and  (12).  These  results  are  given  in  table 
II  for  01=  100  cps.  Four  different  frequencies  were  run,  .01,  1,  100,  and 
1000  cps,  to  determine  the  interval  sums  up  to  t  =  ^ ^ 

equation  (9)  and  t  =  10('^;^)  in  equation  (10).  In  all  cases,  the  sum  con¬ 
tributions  added  up  to  results  accurate  to  four  significant  figures  as 
compared  with  the  final  computed  answers  for  E*  (cu)  and  E”  (oi)  using  the 
full  time  regime.  For  cu=  1  cps,  10(-q}")  =  10  seconds.  Thus  an  accurate 
computation  for  E*  (iou)  could  be  obtained  from  stress  relaxation  tests  run 
for  only  10  seconds. 

The  minimum  time  at  which  the  first  data  point  must  be  known  in  order 

to  compute  E*  (ui)  and  E"  (cu)  accurately  to  four  significant  figures  for 

-3  -5 

cu  =  1  cps  is  1  x  10  sec.  for  E*  (cu)  and  1  x  10  sec,  for  E"  (cu) , 

The  time  range  over  which  the  stress  relaxation  modulus  for  polyiso¬ 
butelene  must  be  known  to  compute  E*  (ic*>)  accurate  to  four  significant 
figures  for  a  given  frequency  is  approximately  5  decades  for  E*  (ou)  and  10 
decades  for  E"  (cu)  .  This  is  true  for  the  frequency  range  from  .01  to  1000 
cps,  which  is  presented  graphically  in  figure  3.  These  limits  apply  only 
to  polyisobutelene  but  those  for  solid  propellant  should  be  similar. 


- 

See  List  of  References 


9 


10 


Figure  2.  E*  (i^)  Polyisobutelene 
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Figure  3.  Time  Regime  in  Which  E(t)  Must  Be  Known  to  Compute 
E*  (iw)  for  Polyisobutelene 


TABLE  I 


STRESS  RELAXATION  DATA  FOR  POLYSIOBUTELENE 


log,Qt  log,QE(t) 

(HRS)  (PSI) 

.-^SO-OOE  +  QZ  .'5-6-5t>aE-+-0-t 
-.14400E+02  .’56460E+01 
- 1 4 2-00  E + 0-2- -  r  5 -6-4 1 0  E +rT  1 - 
-.14000E+02  .56260E+01 
.l?8ejEt02  . 

-.136100E  +  02  .56060;E+01 

-  -—.-1 34-d<3E+02--  J  55-9-]:0EH--at  - 
-.13200E+02  .55710E+01 

— .-13  Olio  E  +-0-2-  -  r  5-5r5-l  eE-+-0'  1 ' 
-.12800E+02  .55310E+01 

- tl2609E-t'02  «  5-!re-1-0Ctm 

-.12400E+02  .54610E+nl 

-  1 2  2-00  E  +-02-  ■  i  1>^T6  OF'+ll  1  ’ 
-.12C00E+02  .53660E+01 

-  1 1806  E+  0-2  -  i  5^1 1 0-£--<^0 1  - 
-.11600E-^02  .  52360E-^01 
■.11400E-<-02  »01410E^01 
-  .  1 12  00E-I-V.2  .50460E-^nl 

—.1 1-^0  E+-0-2  -  41136  OE-KT 1 ' 
-.  10800E-102  .48160E-101 
■  ‘-.-t06-08E+-0-2-s4-6-8-60E-^^01- 
-.10400E-102  .45560E-101 


■.lOOOOE-102 

■-.-98fK10E+-0-l-- 

-.960008-101 

■-.94-CK10E+-o-lr- 

•.92000E-101 

'^.9ooooeTor- 

■.8.8000E-101 

-i86-0tl0E4-0-t- 

■.84000E-101 

'  i  8-20-0  eE-+-at- 

-.80000E-101 


.42860E-ini 

lA-l-^lOE-i-ni 

.40260E-ini 

i  3-806  OE--i-at- 

.37660E-101 

.  3631gE-iCn 

.34910E-101 

-.713560E--fni- 

.  32160E-101 

-.13O86OE-101- 

.29660E-101 


•.76000E-101  .27460E-ini 
'  i  7-400 0 E^ll  1  -  -.-2-6 4 60E-i  0 1- 


•.72000E-101 

'r7-OO0eE-+Ol- 

■.68000E-101 

■TfrfrflOOeKll- 

■.64n0OE-iOl 

'.62000E+0)r 

■.60000E-101 

-i56000E+Oi - 

-.56000E-101 

-rOOOOOeiOl- 

-.52000E-101 

-.80000E-101 


.25460E-101 
-.-24560^-101- 
.23710E-101 
-.23010£'^0t- 
.22460E-ini 
i22010E+01- 
.21660E-101 
.214108+01-- 
.21210E+01 
->"210108+01- 
.20360E+01 
.  20760P+0-1- 


log.Qt  Iog,QE(t) 

(HRS)  (PSI)'^^ 

-.48000E+01  .20660E+01 
- .  46OO0E+-01--.-2061-O8+O1- 
-.44000E+01  .20560E+01 
420008+01—-.  20510E+01 
-.400008+01  .20460E+01 
- .  38000E  +  ei  ■V204-1-OE+O1- 
-.36000E+01  .20360E+01 
-.34000E  +  01-  -.  20 3 -1-OE  +0'!- 
-.320n0E+01  .20160E+01 
30000E+Oi-  -.-200108+01- 
-.28000E+01  .19860E+01 
-.26000E  +  01  .197-1-OE+ni 
-.24000E+01  .19560E+01 
-.22000E+01  .19410E+O1 
-.20000E+01  .19160E+ni 
■— *  180008+01-  y-18960E+n-l 
-.16n00E+01  .16710E+01 
-.140008+01  .i8460E+-01 
-.12000E+01  .18160E+ni 
-.  10000E+01---i-i-7860E+^n] 
-.800nOE-00  .17510E+01 
-;'6noooE*-ocr  .iTliOF+oi 
-.40000E-00  .16660E+ni 
-.20000E-00  .16160E+-O1 
.OOOOOF-99  .15560E+01 
.20000E-00  .14910E+O1 
.40000E-00  .l4210E+ni 
-r600008-00~rl'3-360€+01- 
.80000E-C0  .12410E+01 
.-ieOOOE+01-  -.-il-3-teE-+--01 
.12n00E+01  .lOllOE+01 
-.-140008+01— r86-100E'--00 
.16000E+01  .7O100E-n0 
•"^i8ooo8+or  .0sitrrT8-Tiri 
.20000F+01  .34100E-n0 
■.22i^OOE+OT-;61-OOOT-ni- 
.24000E+01-. 3390CE-O0 
.26C00E+Or-'.83On0F-nO 

NOTE; 

PRINTOUTS  MAY  BE  PUT  INTO  A  MORE;  READABLE! 
FORM.  FOR  EXAMPLE  THE  PRINTOUT  -.45000E 
+  02  BECOMES  —  45.0,  BECAUSE  THE  E  +  02 
OF  THE  PRINTOUT  DESIGNATES  THE  NUMBER  AND 
DIRECTION  OF  THE  DECIMAL  POINT.  NOTE  THE 
FOLLOWING  examples: 

+  0,10600  E  +  02  -  +  10,6 
'+  0.82000  E  +  01  -  +  8,20 
+  0.40000  E  -  OO  -  +  0,40 
+  0,00000  E  -  99  -  ”  0 
+  0.12000  E  -  01  -  +  0.012 
+  0,23000  E  -  02  -  +  0,0023  pm. 
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TABLE  II 


INTERVAL  CONTRIBUTIONS  FOR  E*(io)) 


T 


n 


n  +  1 


rif 


.99999E-45  .44877E-01  .40039E-02  .'■l80^9E-lb  .40039E-02 
•39810E-14  .44877E-01  .63216E-02  .26963E-10  .23177E-02 
.63095E-14  .44877E-01  .99115E-02  .66131E-10  .35898E-02 

.100q_0E-13  .44877E_-01  .  15439E-01  .16146E-09  .53280E-02 

V i  5 84 8' E -'13  ' 4 Ve" llE-oi".  2'4 0 0 1 Y- 6 1 3 9 6 3*9 E -09"  Y  8 5 6 1 9 E -  0 2 
.  251  18E:-13_  .44877E-nl  .37186E-01  .96695E-n9  .13185E-01 
.  3'9'8ibE'-r3  V4'4Y7yE-nY  ■.'57  2Y8E-'01  2'’3''32bE'-6  Y  V20072E-0l' 

.63095E-13  .44877E-01  .87639E-~01  .55941E-08  »30380E-01 

.lOOOOE-12  .44877E-ni  .13362E-00  .13419E-n7  .45983E-01 
.15848E-12  .44877E-nl  .20243E-00  .31798E-07  .68809E-01 
.25118E-12'  .44877E-01  .3030"5E-00  .73637E-07"  .  10062E-Co' 
.39eiOE-12  .44877E-nl  .44768E-00  .16766E-06  .14462E-00 
.  63095E-12  .■44877E-01  .65318E-00'  .37740E-n6  .20549E-Od' 
.innnOE-ll  .44878E-01  .94161E-00  .83978E-n6  .28863E-00 
.  15848E-1  1  .44B8QE-ni  .  13359E+01  . 18140E-n5  .  39409E-00 

.25118E-11  .448S4E-ni  .1849bE+01  .37430E-05  .51395E-00 

398'iOE-il'  .44891F-'^1  .25043E+01  ."75548E-03  .65452E-00 

.63n95E-ll  .44906E-'~'l  .33242E+01  .  14977E-^4  .ai982E-00 

.  lbonOF-10  .449’’5F-o]  .43?]7E  +  01  .28R5‘^E-n4  .99756E-00 

.  15848E-10  .449n9E-^l  .‘=-5080E  +  01  .54344E-n4  .HB/SPE  +  Ol 

.25n8F-in  .45091F-^]  .69017F  +  ni  .10119E-03  .  1  3937E  +  01 

.39810F-1C  .4527eF-nl  .853r3E+01  .18732E-03  .162a5E+01 

.63095'F-in  .45623E-f^l  .  1042  1E+02  .34481E-n3  .18914E  +  01 

.10no3E-n9  .46248E-"1  .12618E+02  .63472E-n3  .21968E+01 
.■15848E-09  .47440E-01  .1519bE+02  .  1  1822E-C2  .25796E  +  01 
.25118E-09  .49654E-01  .18247E+02  .22138E-02  .30491E+01 

.39810C-09  .'^0776E-^:  .2]829F  +  n2  .41224E-n2  .35824E  +  01 

.63r>95F-99  .tl407E-^l  .26915E  +  02  .76311E-02  .41860E  +  01 

.innonF-08  .75372F-91  .30851E+02  .13964E-ni  .4b353E+01 
.  15848F-nB  .lCr93F-'^r)  .36433C  +  02  .25561C-'^1  .  55320E  +  01 
.251188-08  .147708-^0  .4?a81E+n?  .487758-01  .ft4478E+01 

.398108-08  .23:»a3F-00  . ‘'-'^368E+02  .86124E-^1  .74840E-<-01 

.63^95F-oa  .^96108-^^'  .'^97848+02  .182278-00  .888898+01 
.  100008-07  .  708978-'^''  .70-,56E  +  02  .  31287E-00  .10802E  +  02 
.158488-07  .131bb8+01  .833378+02  .609908-00  .13230E+02 

.251188-07  .25216E+0]  .9/a30E+C2  .120278+01  .16493E+02 

.39810E-0y .49181 E+oT . 1 2050E+03  .239648+01  .20678E+02 
.630958-07  .V7399F+01  .146568+03  .48218E+01  .26051E+02 
.lOOOOE-06  .19555E+02  .17943E+03  ,981588+01  .328718+02 
.158488-06  .39736E+02  .220118+03  .201808+02  .40678E+02 
.2‘3il8E-b6‘y8l201F+02^  .26636E+03  .414658  +  02  .462578  +  02 
.398108-06  .161048+03  .301218+0?  .798458+02  .348478+02 
.‘6Yog’5E-b6’".  2ab^YF+’o?  '  .268448+03  .  119288  +  0?-. 327688  +  02 
.lopOOF-05  .319298+03  .979238+02  . 3897 1 8+02- . 1 705 28+03 
.158488-05  .859038+02  . 7 7392 E+02-. 233398+03- . 205308+02 
.251188-05  .308208+03  .206618+03  .222308+03  .129228+03 
.3981  bE'-d'5".y08948  +  b3’  .2  75308+03-.  99262  8  +  02  .686888+02 
.630958-05  .289428+0?  .847838+02  . 80486E+02-. 190528+03 
.  10006E-0'4”.  225  198  +  0?  . 42 794  E +02-.642 35 E  +  02-.  4 198 88+02 
_ .  1_6 8 4 8 8 - 0 4_j.  821918+02  . 1  8 5  1  58  +  03-.  143008  +  0  3  .142368+03 
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TABLE  II 


INTERVAL  CONTRIBUTIONS  FOR  E*(i,w)  (Cont) 


E'  , 

E’’  , 

E’  .-E’  E”  .-E'* 

n 

n  +  1 

n  +1 

n  +  1  n  n+t  n 

.25118E-04' 

. 19306E+03 

. 15458E+03 

.  n087E+n3-.  3  0  57  0E  +  02 

.39810E-04 

. 19306E+03 

. 15458E+03 

.28476E-03  .  28'^76E-03 

.63095E-04 

. 19306E+n3 

. 15458E+03 

.i7761E-n3  .17761E-03 

. lOOOOE-03 

. 19306E+03 

. 15458E+03 

.11078E-03  .11078E-03 

. 15848E-63 

. 19306F+03 

.15459E+03 

.6910CE-04  .69100E-04 

.25118F-03 

.19306E+n3 

.  ]  5459E+03 

.43102E-04  .43102E-04 

•39810E-03 

. 19306E+03 

. 15459E+03 

.23919E-n3  .23919E-07 

.63n95E-03 

. 19306E+03 

. 1 5459E+03 

.14880E-O3  .145a0E-03 

.  lOOOOZ-02 

.  19306E+n3 

.  15459E  +  03 

.8887’4E-n4  .88874E-04 

.  15848E-02 

. 19306E+03 

.  15459E  +  03 

.54172E-04  .54172E-04 

.25118E-02 

.  19306E+03 

. 15459E+03 

.33020E-04  .33020E-04 

.39810E-02 

• 19306E+03 

.  15459E  +  03 

.20127E-04  .20127E-04 

.63095E-02 

.  19306F  +  03 

.  1  5459E+03 

.33693E-04  .33693E-04 

. inoooE-oi 

.  193n6E+n3 

.  15459E  +  03 

.1291-8E-04  .  12918E-04 

. 15848E-01 

.  19306E+03 

.  15459E+03 

.12n93E-n4  .12093E-04 

.25118E-01 

.  19306E  +  O3 

.  1  5459E  +  03 

.72041E-n5  .72041E-03 

.39810E-01 

. 19306E+03 

.  1  5459E+03 

.63445E-05  .61448E-03 

.63095F-01 

.19^06F+03 

ili459£+03 

._36I83E-05  .a61R3E-Q8 

.  lOOOOE-00 

.  19306E  +  03 

.  15459E  +  03 

.28837E-i^8  .28837E-0=> 

.  15848E-00 

.  19306E  +  03 

.  13439E+03 

.21802E-C5  .21302E-05 

.25118E-00 

.  19306F+03 

.  15459E+03 

.15790E-03  .15790E-09 

.39810F-00 

.  193r6F  +  03 

.  15459E+0? 

.110271-03  .11027E-05 

.63095E-0n 

.  19?^6E+''>3 

.15459E+03 

.88313E-06  .83313E-06 

.  inoOOE  +  01 

.  19306E  +  n3 

.  15459E+03 

•  5664dE“*06  •86648E  — 06 

.  18848E  +  01 

.  19306E+03 

.15459E+03 

.3E499E-n6  .334O9E-06 

.2811BE+01 

.193068+0'=^ 

.  15459E+03 

.27631E-n6  .276R1E-06 

.3981OF+01 

.  193n6E+n^ 

.  1  54‘='9E+04 

.  1  7729E-nf,  .17729E-06 

.63O95E+01 

.  19306F  +  0? 

.  15459E  +  03 

.1  1860E-06  .n8f,OE-06 

.9999’9E  +  0V 

M9yo6F  +  03 

.18459E+63 

.68‘43EE-07  .O84'a5E-07 

.15848F+n2 

.  193068  +  03 

.18459E+03 

.4060nc_n7  ./496r'0E-07 

.25118E+02 

.  19306E+0'\ 

.15459E+03 

.24993E-07  .24993E-07 

.39810E+02 

. 19306E+0? 

.15489E+03 

.16393r-n7  .16395E-07 

.63n9  5“E+n2‘' 

yi9306E+03' 

.‘l  5459E'+d3' 

".43889E-0a  .43859E-08 

.99999F+02 

.19306E+03 

.  15459E+03 

.522P3E-08  .82283E-08 

.15848E  +  0'3 

.  193d6E+03 

.1545'9E+63 

■.32316E-08  .R2816E-08 

.25118E+03 

.  19306E+03 

.  15459E+03 

.  .  13373E-ne  .  13375E-08 

NOTE! 

PRINTOUTS  >MlY  re  PUT  INTO  A  MORE  READABLE  FORM. 

FOR  EXAMPLE  THE  PRINTOUT  -  0.045000  E  +  02 

BECAMES  -  45.0,  BECAUSE  THE  E  4-  02  OF  THE 

PRINTOUT  DESIGNATES  THE  NUMBER  AND  DIRECTION 

OF  THE  MCIMAL  POINT.  NOTE  THE  FOLLOWING  EXAMPLES: 

+  0.10600  E  +  02  -  +  10.6 

+  0.82000  E  +  01  ■  +  8.20 

+  0.40000  E  -  00  -  +  0.40 

+  0.00000  E  .  99  • 

0 

+  0.12000  E  -  01  -  +  0,012 

±  0-23000  E  -  02  -  +  0.0023 

ClO-92.3 
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LIST  OF  SYMBOLS 


(jj  ...  circular  frequency 

t  . . •  t ime 

E(t)  ...  stress  relaxation  modulus 
E*(iaj)  .  . .  dynamic  modulus 
E'(cu)  ...  real  part  of  dynamic  modulus 
E’'(cu)  ...  imaginary  part  of  dynamic  modulus 

E^  ...  constant  value  approach  by  stress  relaxation  modulus  for  large 

t 

^  (t)  ...  E(t)-E^ 

m  ...  see  equation  (5) 

a^  ...  see  equation  (10) 

...  sec  equation  (10) 

...  see  equation  (10) 
p  ...  see  equation  (11) 

q  ...  see  equation  (11) 

x\i  ...  see  equation  (7) 

n2  ...  see  equation  (8) 
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APPENDIX  A 


1620  COMPUTER  PROGRAM 

1.  OPERATING  INSTRUCTIONS  AND  INPUT  DATA  FORMAT 


Card,  No. 

1  Tide  card  -  Any  Hollerith  Information  in  first  45 

columns  of  card 


2  Parameter  card 


3  Stress  Relaxation 

4  Table 

5 


-  NT,  NW,  SCL,  FORMAT 

(215, E  14.8) 

NT  =  number  of  time  increments  in  stress 
relaxation  table . 

NW  =  multiple  of  2  ^  at  which  the  pro¬ 
gram  will  begin  using  equation  (12). 

SCL  =  scale  factor  to  multiply  input 
frequency  by  in  order  to  make  the  fre¬ 
quency  and  stress  relaxation  data  in 
compatable  units;  SCL  is  assumed  to  be 
2  7rif  columns  11-25  are  left  blank. 

-  Logj^Q  logj^Q  (E(tj^)) 

Logio  (t2)=  loSlO  (E(t2)) 

Logio  (1:3)  ,  log^Q  (E(t2)) 


M+2 

M+3  Frequencies 

M4-4 

M-l-5 


“  Log^^O  (^m)>  ^oglO 

cjI 

cj2 

(1)3 


M+N42 

M-Hl-i-3 

M-iN+4  Title  card 
M*fN+5  Parameter  card 
etc . 


con 

blank  card  (if  new  stress  relaxation 
table  is  to  be  read  in) 


Format  for  Frequencies  and  Stress 
Relaxation  Table  (2E10 .5) 


( 


\ 


A-1 


2.  PROGRAM  LISTING 
E  STAR  PROGRAM 

DIMENSION  TN(90),  EN(90) ,  AN(90) ,  BN(90) 
TPI  =  6.2831853 
TEN  =  10. 

1  READ  10 

PRINT  10 
PUNCH  10 

10  FORMAT  (45H 

READ  20,  NT,  NW,  SCL 
PUNCH  20,  NT,  NW,  SCL 
20  FORMAT  (2I5,E14.8) 

TRA  =  NW 
TRA  =  TRA*TPI 
IF  (SCL)  28,  24,  28 
24  SCL  =  TPI 

28  DO  30  I  =  1,  NT 

READ  40,  TN  (I),  EN  (I) 

30  PUNCH  50,  TN  (I),  EN  (I) 

40  FORMAT  (6E10.5) 

50  FORMAT  (6E11.5) 

TN  (1)  =  TEN**TN  (1) 

T1  =  TN  (1) 

EN  (1)  =  TEN**EN  (1) 

El  =  EN  (1) 

DO  60  I  =  2,  NT 
TN  (I)  =  TEN**tn(I) 

T2  =  TN  (I) 

EN  (I)  =  TEN**EN  (I) 

E2  =  EN  (I) 

BN  (I)  =  LOGF  (E1/E2)/(T2-T1) 

AN  (I)  =  El*  EXPF  (BN(I)*T1) 

T1  =  T2 
60  El  =  E2 

70  READ  40,  WC 

IF  (WC)  80,  1,  80 
80  PUNCH  50,  WC 
100  W  =  WC*SCL 
110  EP  =  0. 

EPP  =0. 

KSW  ^  0 
TRAW  =  TRA/W 
W2  =  W*W 
T1  =  TN  (1) 

DO  170  I  =  2,  NT 
T2  =  TN  (I) 

IF  (T2-TRAW)  120,  120,  130 


A- 2 


PROGRAM  LISTING  (Cont) 


120  WT  =  W*T2 

CW  =  COSF  (WT) 

CX  =  CW 

SW  =  SINF  (WT) 

SX  =  SW 

125  WT  =  W*T1 

CWT  =  COSF  (WT) 

SWT  =  SINF  (WT) 

TE  =  AN  (I)*W/(BN(I)*BN(I)  +  W2) 
TT  =  EXPF  (BN(I)*T1) 

FN  =  (BN(I)*SWT-fW*CWT)/TT 
FNP  =  (BN(I)*CWT-W*SWT)/TT 
TT  =  EXPF  (BN(I)*T2) 

FN  =  TE*(FN-(BN(I)*SW-fW*CW)/TT) 
FNP  =  TE*(FNP-(BN(I)*CX-W*SX)/TT) 
GO  TO  160 

130  IF  (KSW)  150,  140,  150 
140  CW  =  0. 

SW  =  1. 

CX  =  1. 

SX  =  0. 

KSW  =  1 
GO  TO  125 


150 

BW  =  BN(1)/W 

FN  =  BNW*(EN(I-l)-EN(I))/(l.+Bm-/**2) 

FNP  =  FN 

160 

EP  =  EP+FN 

EPP  =  EPP+FNP 

170 

T1  =  T2 

PRINT  180,  WC,  EP,  FN 

PUNCH  180,  WC,  EP,  FN 

180 

FORMAT  (5H  EP(F8.2,2H)  =  Ell. 5,  12H 

PRINT  190,  WC,  EPP,  FNP 

PUNCH  190,  WC,  EPP,  FNP 

LAST 

TERM  =  Ell. 5 

190 

FORMAT  (5H  EPP  (F8.2,2H)  =  Ell. 5,  12H 

LAST 

TERM  =  Ell. 5 

GO  TO  70 
END 
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DISTRIBUTION  LIST 


Copy  No. 
1 

thru 

10 


11 

thru 

15 


16 


17 


18 


19 


CONVERSION  OF  STRESS  RELAXATION  MODULUS  TO 
DYNAMIC  MODULUS  FOR  A  VISCOELASTIC  MATERIAL 
MTO-270-2 

Recipient '  Mail  S 

Headquarters 
Space  System  Division 
United  States  Air  Force 
Air  Force  Unit  Post  Office 
Los  Angeles  45,  California 
Attention:  SSSC 

Space  Technology  Laboratories 

P.  0.  Blx  1310 

Norton  Air  Force  Base 

San  Bernardino ,  California 

Attention:  Mr.  J.  R.  McGibbeny,  Manager 

Stage  III  MINUTEMAN  Propulsion 
For:  Mr.  R.  W.  Sears 

Space  Technology  Laboratories 

One  Space  Park 

Redondo  Beach ,  California 

Attention:  Mr.  J.  R.  McGibbeny,  Manager 

Stage  III  MINUTEMAN  Propulsion 
For:  Mr.  R.  L.  Greengard 

HP/CMO 

Bacchus  Works 
Magna,  Utah 

Attention:  Lt.  Col.  G.  E.  Scott,  AFOIC 

Chief,  Philadelphia  Contract  Management  District 

United  States  Air  Force 

1411  Walnut  Street 

Philadelphia  2 ,  Pennsylvania 

Attention:  Mr,  S,  E.  Patti 

The  Boeing  Company 
Aero-Space  Division 
P.  0.  Box  3985 
Seattle  24,  Washington 
Attention:  Mr.  T.  A,  Wilson 


Distribution 
Copy  No> 

20 


21 


22 

23 

24 

25 

26 

27 


List  (Cont) 


Mail  Stop 


Mr.  W.  E.  Howell,  Director  of  Contract  Adm. 
in  turn  to:  Mr,  M,  C,  Burgy 
Dr,  L.  G,  Bonner 
Chemical  Propulsion  Division 
Hercules  Powder  Company 
Hercules  Tower 
910  Market  Street 
Wilmington  99,  Delaware 

Mr.  J.  E.  Greer,  Director  of  Air  Force  Programs 

Chemical  Propulsion  Division 

Hercules  Powder  Company 

Hercules  Tower 

910  Market  Street 

Wilmington  99 ,  Delaware 

Mr.  J.  L.  Shrout 

Los  Angeles  Sales  Office 

3460  Wilshire  Blvd. 

Los  Angeles  5,  California 

Mr.  J.  C,  Foster,  Manager 
Allegany  Ballistic  Laboratories 
Cumberland ,  Maryland 

Mr.  N,  L.  McLeod,  Manager 
Hercules  Powder  Company 
Kenvil  Works 
Kenvil,  New  Jersey 

Mr,  B.  H.  Sleight,  Jr.,  Manager 
Hercules  Powder  Company 
Radford  Ordnance  Plant 
Radford,  Virginia 

Dr.  Dwight  F.  Gunder 
463  West  5th  Street 
Loveland ,  Colorado 

Mr.  L.  B.  Johnston,  Manager 
Hercules  Powder  Company 
Rocky  Hill  Plant 
Rocky  Hill,  New  Jersey 


Distribution  List  (Cont) 


Copy  No. 

Recipient 

Mail  Stop 

28 

thru 

J.  M.  Anderson 

117 

42 

43 

W.  M.  Beebe 

117-T-43 

44 

D.  E.  Boynton 

IOO-G-2-1 

45 

R.  N.  Chappell 

117-T-48 

46 

H.  L.  Holt 

117-T-43 

47 

R.  E.  Keeffe 

506 

48 

E.  B«  Magrab 

117-T-43 

49 

J.  H.  Main 

505 

50 

L.  F.  Myers 

117-T-48 

51 

M.  D.  Nay 

117-T-48 

52 

A.  H.  Peterson 

112 

53 

D.  E.  Richardson 

117-T-48 

54 

J.  N.  Sherman 

505 

55 

S.  R,  Swanson 

117-T-48 

56 

T.  L.  Vose 

117-T-43 

57 

and 

Publications  Library 

700 

58 

59  (reproducible) 

Publications  Library 

700 

60 

and 

Library/Central  File 

509-C 

61 

J.  L.  Morse  (letter  only) 

502 

A.  H.  Nielson  (letter  only) 

81-16-2 

J.  L.  Andrew  (letter  only) 

505 

M.  Cohen  (letter  only) 

62 

R.  L.  Griffis  (J.  R.  Larson) 

700 

